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ABSTRACT 


A study is made to determine the optimum shape of an artillery 
shell used to deliver improved ammunition. The "optimum shape" here is 
the one that minimizes the dead mass of the shell for a given payload; 
equivalently the problem is to minimize the surface of a shell with 
specified volume and total drag coefficient. Linear aerodynamics and 
several simplifying assumptions are used in stating the drag coefficient 
dependence on body shape and flight parameters. Restricting the scope 
of the problem to specific families of curves greatly simplifies the 
calculations and yet provides valuable guidance in determining optimum 
shapes; in particular, cone-cylinders and parabolic tangent-ogives are 
studied. Diagrams are presented showing the dependence of surface area, 
for unit volume, on the parameters 6 (nose thickness-ratio) and k (nose 
percentage of projectile length). It is found that optimum shapes do 
not have large length-to-diameter ratios and are nearly independent of 
Mach number in supersonic flight. 



3 




TABLE OF CONTENTS 


(. 

■I 

r. 

Ji' 



Page 


I. 

IT. 


III. 


IV. 


V. 



ABSTRACT . 3 

TABLE OF CONTENTS . 5 

LIST OF ILLUSTRATIONS . 7 

LIST OF SYMBOLS . 9 

INTRODUCTION . 11 

DRAG COEFFICIENT . 12 

A. Assumptions . 12 

B. Values of the Constants . 14 

CONE-CYLINDER . 15 


$ 



A. Optimal Projectile with Prescribed Drag Coefficient.15 


B. Optimal Projectile with Zero Drag . 17 

C. Cone-Cylinder with Prescribed Diameter . 17 

PARABOLIC-NOSE PROJECTILE . 17 

CONCLUDING REMARKS . l8 

REFERENCES . 28 

DISTRIBUTION LIST . 29 


















LIST OF ILLUSTRATIONS 


Figure Page 

1. Diagrams of Cone-Cylinder and Parabolic-Nose . 19 

Tangent-Ogive Projectiles. 

2. Variation of k and A with 6 for Cone-Cylinders . 20 

with C Do = .255 (M w = 2.5) . 

3. Series of Cone-Cylinders of Unit Volume with Drag . 21 

Coefficient C_. = .255 (M ■ 2.5). 

Do “ ' 

4. Variation of A with 5 Near Its Minimum (for . 22 

C Do * .255) and Diagrams of Five Cone-Cylinders 

in This Region. 

5. Variation of C_ with 6 and M. 23 

Do 00 

6. Series of Cone-Cylinders Fulfilling Condition . 24 

3(A)/3k = 0, with C Do * 0. 

7. Variation of Cone-Cylinder Surface for Fixed . 25 

Volume and Fixed Diameter. 

8. Variation of k and A with 6 for Parabolic-Nose ... 26 

Tangent-Ogive Projectiles for Which Eq. (10) is 

Satisfied. 

9. Series of Unit Volume Parabolic-Nose Tangent- . 27 

Ogive Projectiles for Cp Q * .293 (M^ « 2.0). 











ft 


K 

r 

i 

I: 

* 



A 


c 


1 



S 

* 

s 

u « 

V 

x 


LIST OF SYMBOLS 


constant used in evaluating head drag 
(c DH = c i 5 ^ 

skin friction coefficient, Eq. (3) 


base pressure drag coefficient 


friction drag coefficient 

head (nose) pressure drag coefficient 

total drag coefficient of projectile, Eq. (1) 


diameter of projectile. Figure 1 

F (6, k) » 0 -- drag coefficient constraint, 
Eqs. (6) and (10) 

nose percentage of projectile length. Figure 1 
constants in Eq. (2) 


length of projectile 
length of gun barrel 
flight Mach number 


average pressure in gun 

= p U 2 /2 

00 00 9 



= r (x), radial coordinate of point 
surface 

projectile trajectory range 

Reynolds number based on l, p,,,, U*. 
of the undisturbed air 

total surface area of projectile 
- S - (ird 2 /4) 

flight speed of projectile 

volume of projectile 
axial coordinate 


on projectile 


and viscosity 


LIST OF SYMBOLS 


5 

X 


nose thickness-ratio, Figure 1 
Lagrangian multiplier 
density of the undisturbed air 



I. INTRODUCTION 


Artillery systems are currently being considered for the delivery 
of improved ammunition. In this application, the shell serves only as 
a container for sub-ammunition; also a small terminal velocity is desir¬ 
able. The question arises: is there an "optimum" shell for this appli¬ 
cation? 

Optimization of shell design is a much pursued subject. Usually 
a shell is considered "optimum" in an aerodynamic sense of having the 
smallest drag for a given volume or given length. While the considera¬ 
tion of drag is always of importance in ballistics, it is of primary im¬ 
portance only for missiles with high terminal velocity or large range. 

It is proposed here that minimizing the mass of that part of the 
shell which is not effective in terminal ballistics is of primary impor¬ 
tance for improved ammunition. This requires us to reduce the "dead" 
mass of the shell itself, which, as mentioned above, serves only as a 
container. There is the temptation to increase the warhead volume by 
simply lengthening the shell; this, however, will increase the dead mass 
at least at the same rate as the usable warhead volume. While length¬ 
ening the shell is, of course, the only way to increase volume for a 
given diameter weapon, it is not optimizing in the sense discussed 
above. The large i,/d shapes* are also difficult to spin-stabilize, and 
this renders their use even more difficult. It is worth pointing out 
that minimum dead mass not only reduces the fabrication and logistics 
costs, but also reduces the weight of the total artillery delivery sys¬ 
tem, since the size of the recoil mechanism is very sensitive to pro¬ 
jectile weight. 

On practical grounds, it must be acknowledged that there are es¬ 
sentially only three basic diameters at the disposal of the designer, 
namely, those of most of the existing inventory of gun barrels. However, 
further variation of shell diameter could be accomplished by use of sabots. 


*Definitions of symbols are found in LIST OF SYMBOLS, page 9 
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The wall thickness of the shell is expected to be small compared 
to the diameter; then the mass of the shell container is proportional 
to the surface, and the problem is to find the shell shape that has 
the smallest surface for a given volume. It is also necessary to 
specify the drag coefficient; this, of course, rules out the trivial 
shape (sphere). 

If the terminal kinetic energy of the shell is negligible, the 
appropriate constraint is not the drag but the drag coefficient, which 
for a given range R and velocity must be no larger than 

C Do = P L (Rq) _1 

and is independent of the bore area for a full caliber shell. Here p 
is an average pressure in the gun, L the length of the barrel, and 

2 

q-p^U^/2. Given the free stream conditions, the drag depenis only on 
the body shape, as does the surface area and the volume. 

The problem of minimizing the surface area under the constraints 
of a specified drag and volume constitutes a variational problem with 
isoperimetric constraints. The problem generally can only be solved 
with the aid of a computer; for axisymmetric missiles, for example, 
the result will be an "optimum" meridional curve. However, a survey 
of the major dimensions of such an optimum shape may be obtained by 
considering only families of shapes; consequently,, we shall restrict 
ourselves to cone-cylinders and parabolic nose-cylinders (tangent ogive), 
shown in Figure 1. 


II. DRAG COEFFICIENT 

A. Assumptions 

We shall assume that the major portion of the trajectory is super¬ 
sonic. This assumption is not necessarily in conflict with the assump¬ 
tion of small terminal kinetic energy; in a more sophisticated approach 
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one would, of course, properly weight the Mach and Reynolds number 
ranges of the trajectory. Wc shall further restrict ourselves to 
slender bodies. These assumptions are not at all essential for the de¬ 
velopment of the subject, but they allow the use of closed-form ex¬ 
pressions for the pressure drag. 


The drag coefficient Cp o may be considered as the sum of the co¬ 
efficients of various components: i) the head pressure drag coefficient 

C. ., ii) the base pressure drag coefficient C nD , and iii) the friction 
(Hi Ud 

drag coefficient C^p. Thus 


C Do = C DH + C DB + C DF 


CD 


The head pressure drag coefficient in linearized supersonic flow 
theory may be expressed as 

C DH = K lf + K 2 { r (dr/dx) } x=k£ ^rCM^-n 172 ] (2) 

where the constants (K 1 , K.) are determined solely by the value of 
1 / 2 -’ * ^ 

for M ro = 2 '1. It is seen that Cp H is completely independent of 
Mach number if dtr/dx = 0 at the shoulder and for finite dr/dx at the 
shoulder, C nu is only a weak function of the Mach number. 

Uri / 

; 

The base' pressure drag coefficient depends on the Reynolds number 
the Mach nunjber and the geometry. The theoretical prediction of the 
base pressure coefficient is not yet possible. For Reynolds numbers 
greater tlVan 10 6 , reference 1 reports experimental investigations of 

the base/pressure for cone-cylinder models in a free-flight range. 

i 

For thesye large Reynolds numbers, the base pressure appears to be in- 

/ 

dependent of Reynolds number. We shall use the empirical data of 
reference 1; i.e., 

/ 

/ C DB = ' 
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for both cone-cylinder and tangent-ogive shapes. 

The computation of the boundary layer and the friction drag on a 
spinning shell is also impossible without rather stringent simplifica¬ 
tion. For the present exercise, it may be sufficient to use an empiri¬ 
cal formula which has been used in free flight and wind tunnel work to 
estimate the friction drag. Accordingly, the skin friction is computed 
on the basis of the turbulent incompressible skin friction over a flat 
plate having the same area as the wetted area of the projectile ( ex¬ 
clusive of base area). The use of the incompressible skin friction co¬ 
efficient in the above-mentioned application is somewhat suspect, since 
the friction coefficient is known to decrease due to compressibility 
effects. However, the use of the flat plate formula in this case is at 
best a rough approximation. 

If C- is the friction coefficient for the flat plate referenced to 

f 2 
the wetted area, given by the Prandtl-Schlichting formula 

C f = 0.455 (log 1Q Re)' 2,58 - 1700 Re' 1 , (3) 

then the friction coefficient for the projectile is 

C pF = [4/(md 2 )] C f S, 

where the Reynolds number is based on the total projectile length, 

★ 

and S is the surface of the nose plus cylinder. 

B. Values of the Constants 

For the head pressure drag coefficient, the following approximation 
was made: 

C DH * C 1 

A value of c^ = 0.81 was used for cone-cylinders, based on data in 
references 3 and 4; a value of * 1.2, from Eq. (2), was used 

for parabolic-nose shells. 
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The value of * .00310, based on Re » 3 x 10 6 , was obtained from 
Eq. (3). 

The following table shows the Mach number dependence of Cpg:* 

M -- 1.5 2.0 2.5 3.0 

oo 

C DB — .1900 .1460 .1140 .0952 

For both the cone-cylinders and the parabolic-nose shells, the 
following set of nominal values for total drag coefficient was pre¬ 
scribed, based on cone-cylinder measurements reported in reference 1: 

M — 1.5 2.0 2.5 3.0 

00 

C n -- .349 .293 .255 .221 

uo 

III. CONE-CYLINDER 

A. Optimal Projectile with Prescribed Drag Coefficient 
The cone-cylinder (Figure la) has the surface area 
S = Tr(kt) 2 (6/2) [(1 + (6 2 /4)} 1/2 + (6/2) +2k"* - 2] (4) 

and the volume 

V = tt( 6 2 /4) (k*) 3 [k' 1 - 2/3] (5) 

The drag coefficient constraint is now 

F = -C Dq + Cl 6 2 + (2C f /6)[(l + (6 2 /4)} 1/2 ♦ 2k -1 .- 2] + Cp B = 0 (6) 

It is convenient here to eliminate £, from the problem by taking 
A = SV" 2//3 as the quantity to be minimized. The problem is solved, 
according to Lagrange's method of auxiliary multipliers, by determining 
the 6, k, and A values that satisfy the following two equations plus , 
the constraint relation, Eq. (6): 


8 [A ♦ XF]/96 - 0 
3 [A + AF]/9k « 0, 
15 


( 7 ) 






whore X is a Lagrange multiplier. 

The parameter k must be less than unity for a real projectile. 

This requirement and the drag coefficient constraint place limitations 
on the conditions under which solutions can exist (within the framework 
of the assumptions). Thus in Figure 2a, the dotted curve demonstrates 
that no cone-cylinder satisfying the constraint is possible if an un¬ 
reasonably low Cp Q is prescribed for the projectile. The solid curve 

shows the ranges of f> and k in which a solution might be found for the 
prescribed Cp Q (= .255). 

Figure 2b shows the variation of A with 6 for = 2.5. Figures 

3 and 4 present a set of conc-cylinders of unit volume, including the 
optimal case, that fulfill the condition of Eq. (6). In addition. 
Figure 4 contains a magnified picture of the behavior of A near the 
minimum point. The optimum rone-cylinder has a small cylindrical sec¬ 
tion and may not be stabilizable by spin. For practical purposes, how¬ 
ever, the shell labeled #2 will be satisfactory, for the length of the 
cylindrical section is increased by a factor of 1.68 over that for the 
optimum shell, while A is multiplied by only 1.005. 

For the prescribed set of C Dq values, the shape of the optimum 

cone-cylinder was found to be only weakly dependent on Mach number, as 
the values in the following table indicate: 


M 

00 

6 

k 

(A) . 
'min 

1.5 

.412 

.837 

6.14 

2.0 

.393 

.846 

6.20 

2.5 

.376 

.863 

6.25 

3.0 

.356 

.868 

6.31 


Since the dotted curve in Figure 2a indicates an unrealistic value 
of C Qo , it would be instructive to determine what values of are 

attainable under the assumptions made here. Thus, in Figure 5, C^ o is 
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plotted against 6, by Eq. (6), for the limiting value k = 1; the curves 
clearly show the lower limits to the total drag coefficients. 

B. Optimal Projectile with Zero Drag 

When the gasdynamic factors are neglected, i.e., when zero drag is 
assumed, the problem becomes a purely geometrical one. The minimal sur¬ 
face is found by solving the following equations for 6 and k: 

3A/3S * 0 3A/3k = 0 

A series of cone-cylinders of unit volume, including the optimal 
case, satisfying the condition 3A/3k ■ 0, is shown in Figure 6. The 
optimum shape here differs very noticeably from those restricted by the 
condition of fixed drag coefficient. 

The minimum attainable value for A for all bodies is 4.84, namely, 
that of a sphere. 

C. Cone-Cylinder with Prescribed Diameter 


Here we look briefly at an alternative problem: namely, fixed 
diameter and drag coefficient, and variable volume. The results of a 
computation for a sample case are shown in Figure 7, with results for 
the corresponding prescribed volume problem; ' It is interesting to 
note that^thtf'two minima occur for approximately the same shape para- 
_-meters. 

IV. PARABOLIC-NOSE PROJECTILE 

The tangent-ogive, parabolic-nose projectile has the surface area, 

2 

correct through terms of order 6 , 



SI;"# 






vt- 5 ..spy, 

mM 

- « 


S - it (k£) 2 [{k _1 - (l/3)>« ♦ ($ 2 /4) + (6 3 /15)] 
and the volume 


V - ir(« 2 /4) (k£) 3 [k' 1 - 7/15] 


The constraint, correct through terms of order o , is 
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( 8 ) 


( 9 ) 









( 10 ) 


F = -C Do +Cl 6 2 + (4C f /6)[{k _1 - (1/3)} ♦ (6 2 /15)] + C DB = 0 

Figure 8 presents the variation of k and A with 6 for » 1.5, 

2.0, 2.5, and 3.0. In all cases, the minimum A occurs at k > 1, so that 
there is no realizable "optimum" parabolic-nose projectile for the pre¬ 
scribed values of Cp Q . However, the curves clearly indicate the type of 

shell shape that one should employ to obtain a low value of A. The small¬ 
est A values occur at k = 1, and are summarized in the following table: 

M -- 1.5 2.0 2.5 3.0 

OO 

6 (k=l) -- .335 .315 .306 .290 

A -- 6.10 6.22 6.27 6.38 

A set of shell shapes fulfilling the condition of Eq. (10) for 
■ 2.0 is shown in Figure 9. Clearly, the more "efficient" shapes 

(smaller A) have the smaller values of fc/d. 

V. CONCLUDING REMARKS 

In the foregoing study, linear aerodynamics and a base pressure 
coefficient independent of shape have been used in the analysis. Also, 
the particular choice of cone-cylinder and parabolic-nose shell families 
is one of convenience and the graphs are only a guide for the "optimum 
shell." Nevertheless, one can draw tentative conclusions from these 
results; namely i) the shell of minimum surface is a small £/d shell 
(-3-4) for fixed volume and drag coefficient, and ii) the optimum 
shape in supersonic flight is nearly independent of Mach number. 


8»d/(kt) * 2 tanf e 
a CONE-CYLINDER PROJECTILE 


Max-udj'ux ) 2 



• *d/(kU 

b. PARABOLIC TAN6ENT-00IVE PROJECTILE 


Figure 1. Diagrams of Cone-Cylinder and Parabolic-Nose Tangent-Ogive Projectiles 






a. VARIATION OF k~ EQUATION (6) 



Figure 2. Variation of k and A with 6 for Cone-Cylinders 

u |»k r - nee /u _ a e\ 












3 ■ .364,i/d > 3.68 
A • 6.32 



3 >.368,i/d >3.48 
A>6.28 

OPTIMUM 
| >>.376, i/d >3.08 
A >6.25 


3>.382,i/d • 2.78 
A-6.29 












’omift 


1.5 

2.0 

2.5 

3.0 

.250 

.206 

.174 

.155 



CONE CYLINDER 
k>l 


Mtfl.5 

Mqo-2.0 

M^‘2.5 

M<d*3.0 


Figure 5. Variation of C n with 6 and M^. 
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CONE - CYLINDER -- ZERO DRAG 


< 




Figure 6. Series of Cone-Cylinders Fulfilling Condition 
a(A)/ak ■ 0, with Cp o ■ 0* 







CONE-CYLINDER (satisfying Eq.(6)) 



a FIXED DIAMETER 



V/d 


b FIXED VOLUME 

Figure 7. Variation of Cone-Cylinder Surface for 
Fixed Volume and Fixed Diameter. 
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I-.100 A-9.13,1/d-14.2 




8-2 
A-8.07 
i/d-9.64 



8-3 
A -6.49 
J/d-4.24 



8-.315 
A-6.22 
i/d -3.29 



Figure 9. Series of Unit Volume Parabolic-Nose Tangent-Ogive 
Projectiles for C Qo ■ .293 (M^ ■ 2.0). 
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